Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Statistical equilibrium and ion cyclotron 
absorption/emission in strongly magnetized plasmas 

Alexander Y. Potekhin 1 ' 2 and Dong Lai 2 

1 Ioffe Physico-Technical Institute, Politekhnicheskaya 26, 194021 St. Petersburg, Russia; palex@astro.ioffe.ru 

2 Center for Radiophysics and Space Research, Department of Astronomy, Cornell University, Ithaca, NY 14853, USA 

5 February 2008 



ABSTRACT 

We calculate the transition rates between proton Landau levels due to non-radiative 
and radiative Coulomb collisions in an electron-proton plasma with strong magnetic 
field B. Both electron-proton collisions and proton-proton collisions are considered. 
The roles of the first-order cyclotron absorption and second-order free-free absorption 
and scattering in determining the line strength and shape as well as the continuum 
are analysed in detail. We solve the statistical balance equation for the populations of 
proton Landau levels. For temperatures ~ 10 6 -10 7 K, the deviations of the proton 
populations from LTE are appreciable at density p < Q.lB^f g cm -3 , where B\± = 
B/(10 14 G). We present general formulae for the plasma emissivity and absorption 
cocfficents under a wide range of physical conditions. Our results are useful for studying 
the possibility and the conditions of proton/ion cyclotron line formation in the near 
vicinity of highly magnetized neutron stars. 
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1 INTRODUCTION 



Cyclotron lines are a powerful diagnostic tool for magne- 
tized neutron stars. The detection of electron cyclotron 
features at 10-80 keV in th e spectra of a numb er of 



ergy, 



eB 

hio ce = h = 1.16 Bi4 MeV, 



(1) 



binary X-ray pu l sars 



iHeindl et al]|2004 iTerada et al.l l200d for recent observa- 



Triimper et~al] Il978l : see 



tions) provided direct confirmation that these are neutron 
stars endowed with strong magnetic fields B ~ 10 12 - 
10 13 G. Numerous theoretical works have been devoted 
to the physics and modelling of electron cyclotron line 



Wasserman & SalDeterl 


1980 


; iMeszaros & Nagell Il985t 


Burnard. Klein fc Aronsl 


198S 


; lLamb. Wang & Wasserman 



1990; Wang, Wasserman & Lamb 1993(; lArava fc Harding 
19991 : 1 Arava-Gochez fc Hardindl2000l ). 



There has been growing evidence in recent years for the 
existence of neutron stars possessing superstrong magnetic 
fields, B > 10 14 G. In particular, soft gamma-ray repeaters 
(SGRs) and anomalous X-ray pulsars (AXPs) are believed 
to be magnetars, whose radiatio n is powered by the decay of 
superstrong magnetic fields (seelThompson fc Duncanfl995l . 



lies outside the X-ray band, but the ion cyclotron energy, 

huj ci = = 0.635 (Z/A)B 14 keV, (2) 

rmc 

lies in the detectable range for X-ray telescopes such as 
Chandra and XMM-Newton when Bu > 0.4. In Eqs. fTJ 
and ([2]), Bu = _B/(10 14 G), m\ is an ion mass, and Z and 
A are nuclear charge and mass numbers. In the last few 
years, absorption features at E ~ 0.2- 1 keV have been de- 
tected in the spectr a of several therma l ly emitting isolated 
neutron stars (e.g., Haberl et al. 2004; Ivan Kerkwiik et al.l 
|2004 Ivan Kerkwiik fc Kaplanl l2006h . While no definitive 
identifications of the lines have been made, it is likely that 
some of these lines involve proton cyclotron resonances at 
B ;$ 10 14 G. Somewhat surprisingly, the observed quiescent 
emission of AXPs and SGRs does not show any spectral fea- 
ture, in particular the proton cyclot r on line around 1 keV 
(e.g., lJuett et al.l 
l2003l : lTiengo et al 



20021: iPatel et alj 12003! : iKulkarni et al.l 



20051 ) . This absence of lines may be nat- 



1996; Wo ods fc Thomps on 2005). Several radio pulsars with 
inferred surface magnetic fields approaching 10 14 G have 
also been discovered (e.g.. iMcLaughlin et alj|2003 ). In such 
superstrong magnetic field regime, the electron cyclotron en- 



urally explained by the effect of vacuum polarization, which 
tends to reduce the line width significant ly in the atmo- 
sphere (thermal) emission for B > 10 14 G (iHo fc Lail l200l 
|2004 lLai fc Hdl2003l : Ivan Adelsberg fc Lail 120061 ) . 

There has been some evidence for ion cyclotron lines 
during several AXP/SGR outbursts, e.g., the 6.4 keV 
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emiss ion feature in SGR 1900+14 l|Strohmaver fc Ibrahim! 
120001) , the 5 keV absorption feature in SGR 1806-20 
[ Ibrahim. Swank fc Parkef 120031 ) . and the 14 keV emission 



feature (and possibly also ~ 7, 30 keV features ) in the bursts 
of AXP IE 1048-5937 (|Gavriil et all l2002j ). There was 
also a possible detection of a 8.1 keV ab sorption feature in 
AXP 1RXS J1708-4009 (|Rea et al.ll20o3 ; but see lRea et all 
120051 ). It is possible that these absorption/emission features 
are produced by ion cyclotron resonances in the corona or 
lower magnetosphere of magnetars (although one cannot 
exclude the alternative possibility that they are produced 
by electron cyclotron resonances in upper magnetospheres) . 

To be specific, in the following we focus on the electron- 
proton plasma (Z — 1, A = 1.008, spin = i) and the pro- 
ton cyclotron resonance with energy hu cp — 0.630 B14 keV. 
Generalization to other ions is outlined in Sect. [7] 

While the physics of electron cyclotron line transfer has 
been extensively studied in the context of accreting X-ray 
pulsars (see above) , several basic issues regarding proton cy- 
clotron line have not been properly considered. Since the ra- 
diative electron cyclotron decay rate is many orders of mag- 
nitude larger than the collisional deexcitation rate, electron 
cyclotron reson ance always takes the form of scattering (e.g. 
lMeszaro"s] [l992). For protons, the radiative cyclotron decay 
rate is much smaller, so the situation is not at all clear. De- 
pending on the plasma density, temperature and magnetic 
field strength, true proton cyclotron absorption and emission 
are possible. Previous calculations of proton cyc lotron lines 
from magnetized neutron star atmospheres (e.g., Zane et ail 



120011 ; iHo fc Laill200ll , 120031 ; iPotekhin et "alll2004l ) assumed 
local thermal equilibrium (LTE) population of proton Lan- 
dau levels. We will see that this assumption is not always 
valid in the case of a magnetar. 

In this paper we study systematically the rates for 
collision-induced proton cyclotron transitions in a magne- 
tized plasma. Combining these rates with radiative transi- 
tion rates, we then study the statistical equilibrium of pro- 
tons in different Landau levels, and use the non-LTE level 
population to calculate the radiative opacities and emissiv- 
ities for different photon modes. Our results serve as an 
crucial ingredient for determining the possibility and the 
physical conditions of proton/ion cyclotron line formation 
in various plasma environments of highly magnetized neu- 
tron stars. 

In statistical equilibrium, populations of excited Lan- 
dau levels of the ions are determined by rates of sponta- 
neous radiative decay and by rates of transitions caused by 
radiative and non-radiative Coulomb collisions (note that 
'true' cyclotron absorption can be separated from scattering 
by considering second-order Feynman diagrams which in - 
clude Coulomb interaction; cf. iDaughertv fc Ventural ll978). 
Order-of-magnitude estimates of the rates of these processes 
and their consequences for the level populations are given in 
Sect. [2] In subsequent sections we consider these processes 
in more detail. In Sect. Owe write formulae for rates of such 
transitions in a proton-electron plasma. Radiative transi- 
tion rates and cross sections are considered in Sect. [4] In 
Sect. [5] we analyse population of the Landau states of the 
ions. Opacity and emissivity for the ion cyclotron resonance 
are calculated in Sect. [6] 



1.1 Landau Level Basics and Notations 

Motion of the ions and electrons in the plane (xy) per- 
pendicular to the magnetic field B (assumed to be di- 
rected along z) is quantized in Landau levels with excita- 
tion energies E N ,± = mc 2 [Vl + 2bN - 1] (N = 0, 1, . . .), 
where b = huj c /mc 2 is the relativistic magnetic parameter, 
cue = \Ze\B /mc is the cyclotron frequency, and Ze and m 
are the charge and mass of the particle (e is the elementary 
charge). For non-interacting particles, every Landau level is 
degenerate with resp ect to a position of the guiding centre 
in the xy plane (e.g. jLandau fc Lifshitzlll976l ). The number 
of degenerate states (for fixed N and longitudinal velocity 
v z ) is L x Ly \Z\/(2na 2 a ), where 



he \ 1/2 

— I = 2.5656 x lO -11 ^ 
eB 



-1/2 



(3) 



is the magnetic length, and L x and L y are normalization 
lengths. In addition, excited Landau levels of the electrons 
usually can be treated as double spin-degenerate. In con- 
trast, the Landau levels of the ions are not degenerate, but 
split with respect to the spin, because of the anomalous 
magnetic moments of nuclei. 

For the electrons, the relativistic magnetic parameter 
is b e = huj ce / (m e c 2 ) — B/Bq, where Bq = m 2 c 3 /(eh) = 
4.414 x 10 13 G. For density p < 7 x 10 6 (A/Z)B 3/ 4 2 g cm" 3 
and temperature T <C %Lo ce /kB = 1.34 x 10 10 B14 K, vir- 
tually all electrons reside in the ground Landau level. For 
T > 2.7 x 10" 4 B~ 2 (p Z/A) 2 K (where p is the density 
in units of 1 g cm~ 3 ), the electrons are non-degenerate. We 
shall be concerned with this regime in this paper. 

For the ions, the cyclotron energy is 7iu c i = 
0.635 (Z/A) B14 keV, and their relativistic magnetic 
parameter is bi — Tiu> c \/(rrnc 2 ) = b e Z(m e /mi) 2 — 
0.68 x 10 -6 (Z/A 2 ) B14. We consider the situation where 
the ions are non-relativistic: Nbi <C 1, -Ejv,_l = Ntuuid. 
Obviously it is always the case if the ion Landau number N 
is not huge. 

In the following we introduce a number of notations for 
various kinds of transition rates (i.e., number of transitions 
per unit time per occupied initial state) between proton Lan- 
dau levels N and N' and corresponding cross sections. Here 
we list these and other related notations and the basic rele- 
vant equations for easy reference: 



Ct>ce, C^ci, Ucp 



Tie , , Tip 



II N 



cyclotron frequencies for electron, ion 
and proton: Eqs. Q, (J2]| 
magnetic length: Eq. @ 
number densities of electrons, ions and 
protons in the plasma 
number density of protons in the 7V-th 
Landau level (N = 0, 1, 2, • • • ) 
transition rate per particle on the level 
N with initial longitudinal (i.e., along 
B) velocity v z , in a volume V for 
scattering on a fixed Coulomb centre: 
Eq. ifATt 

corresponding cross section: Eqs. (|A7[) . 
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transition rate for protons scattered on 
those electrons which rest on the ground 
Landau level, assuming the relative ini- 
tial longitudinal velocity v z : Eq. (5) 

partial cross section for such scattering, 
normalized to v z : Eqs. (4), (8) 
transition rate for charged particle 1 
scattered on particle 2 in volume V at 
relative longitudinal velocity v z , under 
the condition that initial and final Lan- 
dau numbers of particle i are Ni and 

respectively: Eq. (|A19|I 
analogous to W NlN2 . N ^ N ^(v z ), but for 
collision of identical particles with even 
(+) or odd (-) total spin: Eq. (|A22|I 

corresponding cross section normalized 
to v z : Eqs. (15]), (16]), (jA19]l 
partial cross section of a proton with ini- 
tial Landau number N2 (regardless of 
its final Landau number) with respect 
to scattering of a test proton from Lan- 
dau number N to A": Eq. (18) 
velocity-averaged partial Coulomb tran- 
sition rate for proton scattering on the 
electrons on the ground Landau level: 
Eqs. @, 0, d 

analogous to , but because of col- 

lisions with other protons: Eqs. (115[) . 

GUI 

spontaneous decay rate: Eqs. (S), (1211) . 

di 

total transition rate due to photoab- 
sorption: Eqs. (2D), (3D) 
total transition rate due to stimulated 
emission: Eq. (|30[) 

cross section of a proton with respect to 
photoabsorption (normalized to speed 
of light) for polarization j (j = 1,2), 
photon frequency ui and direction n: 
Eqs. (24), (25), (26), (47) 
partial photoabsorption cross section 
(normalized to speed of light) for basic 
polarization a (a = 0, ±1): Eqs. (25), 
d>), CD 

total and partial cross sections (nor- 
malized to speed of light) for basic po- 
larization a, for transitions caused by 
free-free photoabsorption of a proton- 
electron pair: Eqs. (33), (34), C3I} 
partial cross section a ai NN' (w) for lon- 
gitudinal velocity v z of the absorbing 
particle: Eqs. gD), (SD) 
cross section of a proton with respect 
to scattering of a photon for basic po- 
larization a: Eq. (37). 



2 ORDER-OF-MAGNITUDE ESTIMATE FOR 
PROTON LANDAU LEVEL POPULATIONS 

In this section we present simple estimates of the relative 
population of protons in the ground Landau level (num- 
ber density no) and the first excited level (ni). Other lev- 
els are neglected here for simplicity, and also for simplic- 
ity we assume that transitions stimulated by radiation are 
unimportant. The cyclotron energy of the proton is hu cp = 
h(eB /m v c) ~ 0.63-Bi4 keV. In this section we use without 
proof formulae for the rates of transitions between proton 
Landau levels, deferring their derivation to the following sec- 
tions. 

(i) Coulomb collisions. The collisional cross section 
involving proton Landau transition N — 1 ~ > A = is 
denoted by a^\v z ), where v z is the relative velocity (along 
the z-axis) between the electron and the ion before collision. 
Detailed balance implies o"jp e '(i; z ) = a^\v' z ), where v z and 



are related by m*v z /2 + hu c 



i*v' z /2, and ra, 



m e mi/(m e + m\) ~ m e is the reduced mass. The collisional 
deexcitation rate per proton is 



pC(pe) 
1 10 



n^(v z a w {v z )) =4V27rn 



(4) 



where n e is the electron number density, (• • • ) de- 
notes averaging over the ID Maxwell distribution 
f(v) oc exp(— m»u 2 /2T), and the Coulomb logarithm 



A 



(pe) 



(to be defined later) is of order of unity for the 



plasma parameters we are interested in. In general, A 



(pe) 



depends on parameter /3 P 



hu cv /T = 73.38 B 14 /T 6 , 



where T is the kinetic temperature for particle motion 
along B, and Te = T/(10 6 K) (throughout this paper, we 
suppress Boltzmann constant, implying the conversion 1 
keV = 1.16045 x 10 7 K). At @ p > 1 we have A^ e) ~ 1 
(and A^ e) ~ f3 p 1/2 \ ki/3 p \ for f3 p < 1 (see Sect. 153). The 



collisional excitation rate is rJ n ' pe - ) = Y^ 1 "* 1 exp(— f3 p ). 
The contribution to Landau excitation from proton-proton 
collisions is of similar order and will be neglected in this 
section. 

(ii) Radiative Transitions. The spontaneous decay 
rate of the first Landau level is 



^C(pe) 



1 10 — 



4 e 2 uj% 
3 m„c 3 



(•») 



We neglect here the radiative absorption and stimulated 



emission. 



(iii) Statistical equilibrium. The relative population 
of protons in N — and = 1 is determined by 



m/no=e-^ [l + rf /rf 



(6) 



When Tfo/Yio "C 1, the Boltzmann distribution (i.e., 
LTE) is recovered, n\ = e~ l3p no. In this case, the cyclotron 
absorption and emission are related by the Kirchhoff law. 



In the opposite case Vi /F 10 <<iC 1, 
rfji/r^o, i.e., collisional excitation (0 — 1 
radiative decay (1 — > 0). 

With n e — p/rrip, we see that 



we have ni/no = 
1) is balanced by 



pC(pe) 

1 10 



28.6p Sr7 /2 Ai pe) 



(7) 

The ratio (7) is larger than unity for ordinary neutron star 
atmospheres, but it can become smaller than unity for mag- 
net ars. 
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The above situation should be contrasted with that 
of electrons. The deexcitation rate of an electron from its 
first excited Landau level to the ground state due to colli- 
sions with protons (treated as classical particles) is of the 
order of 4\/27r n p o„e 4 m e /h 3 for /3 e = huj ce /T > 1 [cf. 
Eq. ©]. The spontaneous cyclotron decay rate of electron 
is rf (e) = 4e 2 w 2 e /(3m e c 3 ), and the ratio Y c /V A is about 
a factor (m e /m p ) 7 ^ 2 smaller than that for protons. Thus 
for electrons, radiative deexcitation is always much faster 
than collisional deexcitation, and there is no true electron 
cyclotron absorption, but only scattering, in the magnetic 
fields of ordinary pulsars and magnetars. 



3 NON-RADIATIVE COULOMB COLLISION 
RATES 

Coulomb collision rates of non-degenerate fermions in a 
strong magnet ic field have been studied by many authors. 
IVentural (1 19731 ) derived collision rates for electrons sc a ttered 
by a fixed Coulomb potential. iPavlov fc Yakovlevl (Il976l ) 
presented transition probabilities for collisions of two non- 
relativistic particles, which interact via a screened Coulomb 
potentia l. As a particular case they recovered the result of 
IVentural l | 1973ft . but in a simpler form. Relativistic expres- 
sions for Coulomb collision rates of non -degener a te fer mions 
in a magnetic field we re derived by iLangerl (Il98ll ) and 
IStorev fc Melrose! (|l987ft . However, since we are interested in 
Landau tr ansitions of ions, we may take the non-relativistic 
approach (jPavlov fc Yakovlevll 1976ft . Accordingly, we do not 
consider Coulomb spin-flip transitions which generally are 
weaker by a factor ~ bi compared to the transitions which 
preserve spin. The spin distribution, however, may affect 
statistical equilibrium through exchange effects. 



3.1 Proton-electron collisions 

General formulae for Coulomb collision rates of two dif- 
ferent particles with arbitrary charges are given in Ap- 
pendix [X27T] Here we consider electron collisions with pro- 
tons, assuming that the electrons remain in the ground 
Landau state. This particular case has been pr eviously 
co nsidered by Miller. Salpeter fc Wassermanl l|l987ft . based 
on IPavlov fc Yakovlev ( 19761 ). For a given relative velocity 
(along z) v z between a proton and an electron, the transition 
rate from proton Landau level TV to TV' is 



W nn>( v z) = n e v z a^ e J,(v z ) = 4ttt 1 n e Qm w£% {u±)/u, 

(8) 

where cr^^,(v z ) is the corresponding cross section, to = 
h 3 I (e 4 m e ) — 2.42 x 1CP 17 s is the atomic unit of time, n e is 
the electron number density, and 



w (pe) , n _ r g t/2 4iv'(*/2) d 



(9) 



Here U± = (u =t u') 2 + u 2 , and u = (m*\v z \/K)a m and 
v! = {m r \v' z \/K)ajn are scaled relative velocities along z, 
which satisfy the energy conservation law u' 2 — u 2 + 
2(TV — N')m*/m p , where m* = m e m p /(m e + m p ) is the 
reduced mass. The parameter u s — k s a m , included in 
u±, is the scaled Debye screening wave number (k^ 1 is 



the Debye screening length). For a neutral electron-proton 
plasma at temperature T we have k s = y/Sim e e 2 /T = 
(1.584 x 10 8 cm _1 )po /2 T 6 " 1/2 . Equations © and © follow 
from (|AT9| and igjOjl of Appendix \K\ with Z x = Z 2 = 1 
and w£x,(u±) = w 0tN: ,n> (n±). Laguerre function I NN i is 
defined by Eq. (|A4| |. 

If the distributions of z-velocities of electrons and pro- 
tons are Maxwellian with temperatures T e and T p , respec- 
tively, which do not depend on the Landau number TV, then 
the relative velocities v z — hk/m* have Maxwellian distri- 
bution 



2ttT 



exp 



2T 



where 



T — (m e + m p )/(m e T p 1 + m v TZ x ). 



(10) 



(11) 



In order to simplify formulae, hereafter we assume 
T e — T p — T, unless the opposite is explicitly stated. 
Then the velocity-averaged partial Coulomb transition rate 

(P e ) 
' JVJV' 



-.C(pe) 



where 

a(p e ) 



4(e 4 /h 2 ) V2 7 rm»/Tn e o 2 n A^ / 
A\f2~n ( m t m p 



A' 



JVJV' 



(pe) 



A' 



Here 



Pv 1V NN' 

du _ 



1/2 



i (P e ) 



(u )g(u)g(u') 



(12a) 
(12b) 

(13a) 
(13b) 



) [with u = u 2 + 2(TV - N')m t /m p ] is the 
step function that ensures the energy conservation, /3* = 
heB /m,cT — /3 p m p /m*, and g(u)g(u') is the correction 
factor, which approximately allows for violation of Born ap- 
proximation as discussed in Appendix IA1. 41 The latter fac- 
tor appreciably differs from 1 only at u < jg 1 ^ 2 , where 
7_b = h 3 B /(m 2 ce 3 ). In the case of electron- proton collisions 
7 " 1/2 w af/s/bl = 0.004848 B~ 4 1/2 . The smallness of 7~ 1/2 
ensures that the approximations used to derive Eq. (|13bf) 



are sufficiently accurate; in this case the 7_g-dependence in 
Eq. (I13b|) is weak (logarithmic). 

By changing integration variable u — -> v! in Eq. (|13b|) . 
and taking into account that w^^j,(u±) = i/j^f^ (tt±), we 
can check that A ( ^ e J, = e ^-N') Afet , and thus 



"-jvjv 



^ p (N-N') r C(pe) 
i JV'JV • 



(14) 



*■ JVJV' 

FigureQ]depicts the function A^!?, (TV > TV') for transi- 



tions between different low-lying proton Landau states, cal- 
culated assuming u B — * and 7^ — ► 0. Figure [2] presents 
the same functions for u s — 0.5. Note that u s = fc s a m = 



4.1 x 10 



-3 1/2 T - 



1/2 R -l/2 



So u s 



0.5 corresponds to 



a rather high plasma density. We see that at any B and 
u s transitions between neighbouring states (TV — TV' = 1) 
strongly dominate. 

Representation (|12b[l is most convenient when j3 p ^> 1, 
because in this case the exponential function in Eq. (I13b|l 
varies much faster than u' and 



,,(p e ) 



and it can be inte- 



grated separately. Hence A 



(pe) 



JVJV' ' 

approaches a constant. 
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In the opposite limit /3 P <C 1 it may be more conve- 
nient to use Eq. (|12a|l , because in this case A^^, is a slowly 
varying function of j3 p . At the first glance it may seem un- 
physical that Eq. (|12ap contains factor which goes to 
infinity as B goes to zero. However, it has a simple explana- 
tion. As long as the Landau numbers of the electron (equal 
to zero) and proton (TV, TV') are kept fixed, Eq. (|12a|l de- 
scribes the partial rate of the collisions in which the transfer 
of the kinetic energy of the motion transverse to the field, 
TV — N'\huj C p, decreases linearly with decreasing B. In the 
classical picture this corresponds to collisions with impact 
parameters increasing oc a m , for which the cross section in- 
creases according to the Rutherford formula. The divergence 
of the classical cross section at large impact parameters is 
eliminated if one takes into account the screening of the 
Coulomb potential. It is also the case for the quantum cross 
section. Indeed, u± in the denominator of Eq. ([9| in general 
includes the term u 2 — feja^- Thus w~r% («±) (and hence 



(pe) 



becomes oc (fc s a m ) 4 when the magnetic length a n 



A 

is much larger than the screening length k~ 



3.2 Proton-proton collisions 

Now let us consider proton-proton collisions. This case is 
more complicated than the previous one in two respects: 
first, there is the exchange interaction described in Ap- 
pendix IA2.2I and secondly, both colliding particles can 
change their Landau numbers (neither of them is confined 
to the ground state). 

Let tin be the number density of protons in the TVth 
Landau state, and let and fj^ be the fraction of such 
protons with spin along and opposite to the field direction, 
respectively (/^ = 1). Then the average rate of proton 
transitions from level N to level N' due to the Coulomb 



collisions is 



r C(pp) _ 
L NN' 



N 2 N' 2 

HfUk + fhfk ) . (is) 



where ( v zO'^ 1 ^ 2 ± N , N i ) * s ^ ne probability, per unit time, that 
two protons in unit volume, which have initial Landau num- 
bers JV and N2, make a transition to the state where they 
have Landau numbers N' and N2 , under the condition that 
their spin projections to B are the same (sign — ) or opposite 
(sign +). The factor | at the sum allows for the quantum 
statistics of identical particles. 

For Maxwell distribution (|10p with m* = m p /2, using 
the results of Appendix I A2I we obtain 

f°° du _ „2 r , . , . 

Xy '[i/' 6 [ W NN 2 ;N>N! ! {U + ) +W NN2 . n , n ^{U-) 

± 2w NN 2 ;N>N 2 ( u -, u +)], (16) 

where u% = (u ± u') 2 + u 2 and u' = [u 2 + N' - N + N 2 - 
N 2 ] 1/2 . Functions w NN2 . N , N , 2 and w x NN2 . N , N ^ are given by 
equations (|A20[l and (|A24[I . respectively. 

Equations (|15p and (|I6[) can be written in the form 
analogous to Eq. (|12b|) . 



(17) 



N 2 



to m e L 



A (pp) 



NN 2 ;N' 



(/jV /iv) (/iV 2 fk ) ^JVAr 2 ;'iV'] ' 



(18) 
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0.01 0.1 1 10 



Figure 3. The dimensionless function N , (Eq. I19a0 with- 

out screening (u B = 0) for the transitions between the proton 
Landau states marked near the curves. Initial Landau number of 
the second proton equals 2Va = (solid lines), JVjj = 1 (dashed 
lines), or N2 = 2 (dot-dashed lines). 



where 

A&SU' = V^E r ^^ u2 e{u' 2 )g(u)g(u') 

x[wnn 2 :N>n^(u+) +u>NN 2 ;N'N! 1 (u-)], (19a) 
A^J N , = 2 ^ ^ e"^" 2 S^' 2 ) g(u) 9 (u') 

N' J ° 

"2 

xw x NN2 . N/N ^{u^,u+). (19b) 

In these equations u' and u± depend on N^- The factors 
g(u)g(u'), with <?(w) defined by Eq. (|A16[) . account for the 
correction due to violation of Born approximation, as dis- 
cussed in Appendix IA1.4I However, unlike Sect. 13.11 here 
7 fl 1//2 = afm p /(2m e v / 6e) = 4.45 B^^ 2 is larger than 1 
for B < 2 x 10 15 G, which reflects the fact that the pro- 
tons are moving much slower than the electrons, therefore 
Born and adiabatic approximations (see Appendix lA1.4[) are 
less applicable to the proton-proton collisions. Nevertheless, 
we use these approximations, considering them as order-of- 
magnitude estimates, which is justified because the whole 
effect of the proton-proton collisions on statistical equilib- 
rium is not very significant, as we will see below. 

In Figure [3] we show Aj%% . N i for the case of negligible 
screening (u s = 0) and 73 = 1, for a few initial and final 
proton Landau numbers N, N', and for initial Landau num- 
bers of the second proton N2 = 0, 1, 2. The decrease of the 
displayed functions at f3 p 3> 1 results from the correction be- 
yond Born approximation. This indicates that an accurate 
evaluation of the proton-proton collision rates would require 
non-Born quantum calculations, which are beyond the scope 
of the present paper. 

In figure [4] we compare some of the curves from fig- 
ure [3] (solid lines) with the case of non-negligible screen- 
ing, fc a a m — Us ~ 0.5 (dot-dashed lines), which can be rel- 




*" cp /T 

Figure 4. The functions N , (solid lines without screen- 

ing and dot-dashed lines with screening parameter u B = 0.5) and 
^iVWV N' ~ ^Nn'-N' (^ asnet ^ lines , without screening), repre- 
senting the expression in the square brackets of Eq. H18H for ran- 
dom spin orientations (/T = f*> = 0.5) and for aligned spins 
(/T = 1), respectively, for transitions N — » N — 1 and N — > 0, 
with N = 1, 2, 3, 4, N2 = 0, as functions of f) p = huj cp /T p . 

evant at rather high density. Also in this figure we show 
a comparison of A]^^.^, (solid lines) with the difference 

A Ni\r -JV' ~ ^nn^n' (dashed lines) which enters Eq. 
when the proton spins are all aligned in the same direction. 

Let us note that if T e 7^ T p , then in this section T p 
should substitute T defined by Eq. (|lip . In particular, (3 P = 
Tiuj cp /T p in Eqs. (JTSJ and JT9J. 



4 RADIATIVE TRANSITIONS 

4.1 Radiative transition rates in magnetized 
plasmas 

Magnetized plasma is a birefringent medium. Electromag- 
netic radiation propagates through it in the form of two 
normal polarizatio n modes j = 1, 2 with polarization vec- 
tors e J (uJ,n) (e.g.. lGinzburgiri970l ). Here, cu is the angular 
frequency and n the unit vector along the wave vector. Con- 
sequently, radiative transition rates depend not only on u, 
but also on j and n. 

Let T^ N> , T^ N i , and F^n' be the rates of transitions 
from level N to N' due to spontaneous emission, photoab- 
sorption, and stimulated emission, respectively^ These rates 
are the total (for both polarizations, integrated over angles 
and frequencies) transition probabilities per unit time for 
one occupied initial quantum state. They can be expressed 
through Einstein coefficients ^jvjv' andBjvjv' (for emission) , 

1 Generally, N and N' may take any values. For instance, free- 
free photoabsorption is allowed for N' $C N (i.e., the photon is 
absorbed while the proton makes a downward transition) as well 
as for N' > N. 
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or B NN i (for absorption). These coef ficients have different 
defin it ions in the liter ature (e.g., cf. Rvbicki fc Lightmanl 
1 19791 . iGinzburd Il97(j and IZheleznvakovl Il996l ). We de- 
fine A NN i, B NN i, and B NN i from the conditions that the 
number of quanta with angular frequencies in the inter- 
val duj and wave vectors in solid angle element dn spon- 
taneously emitted by a unit volume during unit time equals 
unA nn i dwdn, and the number of quanta emitted or ab- 
sorbed under the action of radiation with the specific inten- 
sity 7^ equals tinBnn'Iiu d a;dn or unB nn i da;dn, respec- 
ti velv |Zheleznvakovlll996r ). This defini tion (or a similar one 
m ]Gmzburdll97oT but not the one in iRvbicki fc Lightmanl 
1 19791 ) is relevant in a strong magnetic field, where the emis- 
sion is neither isotropic, nor unpolarized. Then 



i ww 



i w/v' 



/ dn / dwB WA r'(j>,n)/ u ,j(n), (20) 
/ dn / duA NN ,{j,cv,n), (21) 

3=1,2-' J 



and the expression for Tff N i is completely analogous to 
Eq. |[2Q]|. 

The quantities A NN i, B NN i and B N / N are related by 
the Einstein re lations (which i nclude polarization depen- 
dence, see, e.g.. lGinzbur3ll970l ; Izheleznvakovll 19961 ) Fl 

Bnn' =B n / n , A nn i = B N i N . (22) 

From the first of these relations, it follows that the stimu- 
lated emission rate is equal to that of photoabsorption with 
interchange of the initial and final levels: Tff N i = ^n'n- 



In Eq. (1221) we have neglected the difference of the group 
and phase velocities of radiation. Einstein relations with al - 
lowa nce for this differenc e are given, e.g., bv lGinzburgl(|l970r ) 
and lZheleznvakovl l|l996h . Note, however, that this difference 
would lead to appearance of additional factors not only in 
Eq. (|22[) . but also in the expressions for photoabsorption 
cross sections discussed in Sect. 14.21 below. 

Spontaneo us cyclotron decay rates have been 
derived by [Paughertv fc Ventura! lfl977ri fsee also 



Daughertv fc Venturalll978l; iMelrose fc Zheleznvakovl | l98ll ; 
Pavlov et alj Il99ll ; iBaring. Gonthier fc Harding! 120051 . and 
references therein). In the non-relativistic limit (Nbi <C 1), 
the decay rates are proportional to bf~ N +1 , multiplied by 
a combinatorial factor. Although the latter factor may be 
large for large N' - N, transitions N — > N' = N - 1 still 
dominate in the non-relativistic regime. The rates of the 
latter transitions are T n 



,22 

4 e u> cv 



3 m v c s 



(23) 



is the natural width of the proton cyclotron line. 

Spin-flip transitions for protons are unimportant in the 
non-relativistic limit, because their rates contain an addi- 
tional factor h comp ared to the dominant transiti ons pre- 
serving spin (e.g., cf. I Melrose fc Zheleznvakovlll98ll ). 

2 Einstein relations i ]22[ i come from the detailed balance 



equation n N i 



n NA NN / 



2 B " 



complete thermodynamic equilibrium, where 
(huj 3 /8tt 3 c 2 ) (e fi "/ T — and the requirement that the coeffi- 

cients A and B must be independent of T. 



4.2 Relation between emission rates and 
photoabsorption cross sections 

The Einstein absorption coefficient B NN i is given by the 
relation 



-Bivjv (j, n) = a NN , (j, w, n)/hu>, 



(24) 



where a NN i is the partial photoabsorption cross section re- 
sponsible for the N — ► N' transition. Equation (|24|) directly 
follows from the definition of B NN i in Sect. 14.11 Together 
with Einstein relations, it allows one to express spontaneous 
decay rates r^ N , through partial photoabsorption cross sec- 
tions. 

In the 'rotating coordinates' (e.g. Meszaros 1992]), the 
polarization vectors e 3 (ui,n) of two polarization modes j = 
1,2 have the components e> a , a = 0, ±1. In the dipole ap- 
proxima tion, photoabsorption cross sectio ns can be written 
as (e.g.. IVentura. Nagel fc Me szaros 19791) 



cr(j,u,n) = ^ a a (uj) |e^(u;,n)| 5 



(25) 



where the component e_ is responsible for the electron cy- 
clotron resonance, and e+ for the ion cyclotron resonance. 
Using equations (|21[l . (|22[) and (|24l) . one obtains 



2 duo 

— a NN>{j,u,n). (26) 



i=i,2 

As already stated in Sect. 14. ll we neglect the difference 
of the group and phase velocities of radiation . This is equiva - 
lent to the 'semi-transverse approximation' (|Venturalll979l ). 
where refraction indices are close to 1, and n- e? « 0. In this 
approximation, the relation Y2j=i 2 ~ 1 holds, where 



A 3 = — 

A a - ^ 

Then equations (|25f) and l)26p gi 



^— j |e^(ix),n)| 2 dn. 



- N'N 



E 



ui dus . . 

o22 (J a,NN l (W). 

6-k z cr 



(27) 



(28) 



Let us suppose that transition N — + N' corresponds 
to an absorption line with a profile (j>(uj) (J 4>{ijj)Aijj = 1), 
which at ui = u)o has a sharp peak with a characteristic 
half- width v <C u^o for the polarization a. At uj ~ c^o, 
let us write the photoabsorption cross section in the form 
a a,NN'( UJ ) = 2^o- aiJV ]v'0(^). Then Eq. ((2"ci)) gives the spon- 
taneous emission rate 



2va n 



NN' <^0 



Let us assume, in addition, that I Ui j(n) can be re- 
placed by the average over the angles under the integral 
in Eq. (|20|l . which we denote Iu,j (this replacement is ex- 
act in the diffusion approximation). Then, using equations 
((24)) and (f27)) . we obtain the photoabsorption and stimulated 
emission rates 

1 r ' c2 r^,v E Mi Jj-. (30) 



3=1,2 



where 



: <f>j(uj) du>. 



(31) 
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4.3 Cross sections at the proton cyclotron 
resonance 

Although cyclotron emission rates can be calculated in the 
framework of the first-order perturbation theory, this the- 
ory is not suitable for the determination of the frequency 
dependence of photoabsorption cross sections and opaci- 
ties. The reasons for that, and the conditions where the 
first-order process still can b e important, were discussed by 
iDaughertv fc Ventural |l978), who stressed that the spec- 
tral dependence of the absorption coefficient is properly de- 
scribed by second- and higher-order processes. Indeed, be- 
cause of kinematic requirements (energy and momentum 
conservation), the first-order absorption is possible only at 
a single frequency at any given angle of incidence. Thus one 
must take into account level broadening in order to obtain 
the spectral absorption coefficient. The broadening is caused 
by the finite life time of the proton in a final state after an 
absorption event. In a macroscopically homogeneous plasma 
this life time is limited only by (1) spontaneous emission 
and (2) interactions with other particles. It is the emitted 
photon in the first case or another plasma particle in the 
second case that carries the energy and momentum needed 
to restore the kinematic balance. Thus a quantum descrip- 
tion of the absorption line shape requires at least two- vertex 
Feynman diagrams. The second vertex may correspond to 
the emission of the photon [case (1); cf. Figs. 1 and 2 of 
IDaughertv fc Ventura|[T978l ] or to Coulomb interaction with 
a charged particle [case (2)]. The first case is scattering, and 
the second is free-free photoabsorption. We now consider the 
cross sections for these two processes for the polarization 
component a = +1, corresponding to the proton cyclotron 
resonance. 



4-3.1 Free-free absorption 

The free-free photoabsorption cross section at any frequen- 
cies and polarizations is given by Eqs. (|B1[) - (IB3|I of Ap- 
pendix|B] It can be presented as a sum of terms correspond- 
ing to transitions of a proton from level N to N': 

= ^ZfN^Z^,NN'(^), (32) 
N N' 

where is the fraction of protons in Landau state TV, and 

ff i \ _ 47re2 u2 E„/SE„'"I;L' 1 jv'H 

c a ,NN'(u) — 7 ; T57 r^— 5-^7 — r, 

m e c (uj + auce) (w — au) cp y + u> z u£(uj) 

(33) 

where f„ is the fraction of the electrons in Landau state n, 
and u a {u>) is a damping factor. The separation of erf into 
a a nn' 1 expressed by Eq. (|32|l , is useful for calculation of 
non-LTE emissivity (see Sect. [6}. 

At uj ~ uj cp and a = +1, there is a resonance: 



with v+(ui) given by Eq. (JB2J . From Eq. 



taking into 



47re 



m p c (ui — LJc P ) 2 + u 2 



(34) 



where u = (m e /m p )i> + (uj cp ) = u p + u p , with u p = u p (ui cp ) 
the radiative damping rate r r!P /2 (note that it could also 
include other damping mechanisms not related to electron- 
proton collisions, such as the damping rate due to collisions 
with neutral particles), and u p the damping rate due to 
electron-proton collisions: 

Up = (m e /m p ) v + (u C p), (35) 



account the condition u <C ui cp , we obtain the rate of N' 
N transitions caused by the resonant free- free emission, 



r ff ~ 4 

1 +,N'N ~ — 



1/2 



«.- 



V A 

STY J (uj — UJ cp ) 2 + v' 1 



duj, 



(36) 



where A 



ff,- 

n,N: n' ,N 



, is given by Eq. (|B10|) . This result is written 
m the form similar to Eq. (I12b|) for easy comparison, which 
shows that by order of magnitude r ff /T c(pe) ~ Y I>p /v. The 
damping factor v is discussed in Sect. 14.3.21 here we note 
only that the ratio r ff /r C( - pe - ) cannot be large, because v ^ 

r,,,/2. 



4.3.2 Scattering 



The resonant cyclotro n scattering (jCanuto. Lodenquai fc Rudermanl 
Il97ll : IVentural H979I ) is a second-order process which is 
common for electrons in white dwarfs and magnetic neu- 
tron stars, and for ions in magnetars. The photon-proton 
scattering cross section is 



t7 + = OTp 



(lj — UJc P ) 2 + V 2 



(37) 



where <tt p = S7ve 4 /(3m 2 c 4 ) is the Thomson cross section for 
protons. 

The determinat ion of the effective damping factor v 
(not considered by ICanuto et al.l Il97ll ) is not trivial. In 
general, this task requires a non-perturbative tr eatment 
|Cohen-Tannoudii. Dupont-Roc fc Grvnberd Il998l) . which 
goes beyond the scope of our paper. However, u(u) can be 
found from the correspondence to the classical physics. 

The naive estimate of v as the sum of total half-widths 
of two Landau levels would lead to replacement of Eq. (|37[1 
by a sum of different Lorentz profiles for different proton 
states N. However, this estimate is incorrect, because it ig- 
nores the cohere nce of equally spaced quantum states, as dis- 
cussed, e.g., by ICohen-Tannoudii et all (| 1998f ) for the case 
of interaction of electromagnetic field with a quantum oscil- 
lator. Interference of transition amplitudes between different 
states leads to the common damping factor (which proves 
to be equal to the classical oscillator damping factor) for all 
transitions which have the same resonant frequency. Thus 
we should put in Eq. (|37[1 the same damping factor as in 



Eq. J3SJI, 



u p + u p at lu 



The frequency depen- 



dence of u is su ggested by analogy with a classical oscillator 
l|Jacksonlll975l ): 



with 



u(lj) = u p (uj) + (m e /m p ) u + (u), 



u p (u}) 



3 m p c 3 



(38) 



(39) 



Thus we recover the dam ping factor that wa s previously 
given without discussion bv lPavlov et all (119951 ). Obviously, 
at the resonance, u p (ui cp ) — u p = T t , p /2. 

Note that for damping of free-free photoabsorption 
(Eq. [B3| l we should include, beside u p (u>), also u e {w) = 
(2/3)(e 2 /m e c 3 ) uj 2 . Then the terms containing factor a in 
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Eq. (|B3P cancel out, and it simplifies to 
2 



3 m*c 3 



LU 2 + V%(U). 



(40) 



5 STATISTICAL EQUILIBRIUM OF PROTON 
LANDAU LEVELS 

5.1 Two- level system 

The model in which only two quantum levels participate in 
the radiative and collisional transitions is helpful for under- 
standing the main features of line formation and transition 
rates. This simplest model can be applicable to the forma- 
tion of the proton cyclotron line if the ground Landau level 
is much more populated than excited ones. For electron cy- 
clotron lines, a similar model was considered previously by 
iNaeel fc Ventural (|l983l 1. 

The statistical equilibrium of two proton Landau levels 
is given by the equation 



(41) 



Let us first consider the case where proton-proton collisions 
are unimportant. Then, taking into account Eq. (|30|l and 
the relation r^ (pe) = rf (pe) e- /3 " (Sect. GEE), Eq. gj) can 
be written in the form 

ni =e -p p _l + efl/(l-e-*) 
no 



l+R + eR/(ePr - 1)' 



R- 



(42) 



where the parameter e = 2£] . A 3 + Jj/B u (at u = u) c 



char- 



acterizes the ratio of the effective radiative energy density 
in the line to its equilibrium value. 

Let us mention three important limiting cases. 

• When either R <C 1 or e = 1, the Boltzmann ra- 
tio ni/no = e _/3p is recovered. This is the LTE situation, 
where absorption and emission coefficients are related by 
the Kirchhoff law. 

• In another limiting case, where R S> 1 and eR <C 1, 
ni/no = e'Pp/R = roi/IYp, that is, excitation of the level 
N = 1 is collisional, but its deexcitation is radiative. This is 
the case for which emission is most prominent. 

• In the third limit, where e <C 1 and eR 3> 1, the 
level N = 1 is excited by absorption of radiation and 
deexcited by spontaneous emission. In this case m/no = 
(8ir 3 c 2 /huj 3 ) J2 A j + J h and 
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Figure 5. Population of the proton Landau level N = 1 relative 
to N = as function of mass density for different values of B and 
T. 



(N 2 ) 



3»V<«.*jiw)}/(Ixb + r S i +rf (pe) ), the factor 
1} is given by Eq. (|18[) . an d Xi is the solution (|42[) . 



(0,1) 



0. 



where c 

which is reproduced when c 

In Figure [5] we show the relative populations of the ex- 
cited proton Landau state (N = 1) as function of density, 
according to Eq. <|43f) . for several B and T values. Here we 
assumed that stimulated transitions are unimportant and 
set e = in Eq. (|42|l . At high density, all curves tend to 
their LTE limit ni/no = e _/3p . Radiative decay rates dom- 
inate at lower densities, where the excited level becomes 
depopulated. 

For the considered plasma parameters, the rates of tran- 
sitions due to proton-proton collisions are of the same order 
of magnitude or smaller than those due to electron-proton 
collisions. Therefore a neglect of the pp rates does not sig- 
nificantly affect the statistical equilibrium. For instance, in 
Fig. [5] this neglect would change ni/no by less than 6 per 
cent. Thus statistical equilibrium can be approximately eval- 
uated with only proton-electron interactions taken into ac- 
count. 



^ J diMip (J, U!, n) = np ^ j dnBoi(j,^, n)/ w j. 5.2 Multilevel system 



Then the spectral power of spontaneous emission is identical 
to that of absorption, and both processes can be treate d 
as non-coherent scattering l|Venturalll979l ; lMeszarod[l99l) . 
Such situation is most common for the electron cyclotron 
absor ption and emission in st rong magnetic fields of neutron 
stars (|Naeel fc Venturalll983l ). but it is not so usual for ion 
(proton) cyclotron processes. 

Taking into account proton-proton collisions, from 
Eq. (|4l"Tl we obtain 



The statistical equilibrium of proton distribution over Lan- 
dau levels is determined by the balance of the total rates of 
transitions from and to every level N, 



n N 



+ £01 



N'<N N'j^N 

+ n N> (r 



£ n N > r 

N>>N 



A 

N'N 



N'N T 1- N'NIi 



(44) 



Here 



21 = I ([( 

no 9-l,,WU 1 



1 4- r (0) 
1 + c 10 



- Xl) 



n ° 2 + 2ci u 
+4 (1 + c«) (xi + eg?)] 1/2 - (1 + $ - 4V - x x )}, (43) 



supplemented with the condition X]jv n iv = 
r^jv = T^jv' + ^nn^ ■ This system is non-linear, because 
according to Eq. (|15p r^^' depends on the distribution of 
Wat. We solve Eq. (|44[) iteratively, starting from the Boltz- 
mann distribution of njy. 
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Figure 6. Populations of the levels N = 1,2,3,4 relative to N = 
as functions of mass density for T = 10 7 K and two values of 
B. The solid lines show the population of the N = 1 level, dot- 
dashed lines N = 2, long-dashed lines N = 3, and short-dashed 
lines N = 4. The dotted lines show n±/no based on the two-level 
approximation. 



An example of the numerical solution of Eq. (|44p is 
shown in Fig. [5] As in Fig. \S\ here we have neglected F B 
relative to F c . The curves of different style correspond to 
the values of njv/no as functions of p for different N from 
1 to 4. The dots correspond to n\/no in the two- level ap- 
proximation (cf. Fig. [5} ■ We see that they coincide with the 
multilevel solution for n\/no (solid lines) within graphical 
accuracy. 
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Figure 7. Upper panel: The opacity components for the polar- 
ization a = +1, as functions of the photon energy, for p = 0.01 
g cm" 3 , T = 10 7 K, and B = 10 14 and 10 15 G (as marked near 
the curves). Dashed lines — ion scattering opacity, dotted lines 
— electron scattering opacity, dot-dashed lines - free-free absorp- 
tion contribution, solid lines - the total. Lower panel: relative 
emissivity, Eq. J53}, for the same plasma parameters. 



6 OPACITY AND EMISSIVITY 

6.1 Relation between emission and absorption 
coefficients 

For each polarization component, the photoabsorption coef- 
ficient can be presented in the form 

= Vnn>(u), (45) 

JV,JV' 

p NN ,(cu) — n N J ^(vz) o- NN ,(v z , u)dv z 



n N > \ ^ r N'('vi)o'NN'{vz,u!)dv' z , 



(46) 



where a NN i (v z , uS) is the (free-free) partial photoabsorption 
cross section for ions in the Landau state N having lon- 
gitudinal velocity v z and going to the final state N' , and 
J~n(Vz) is the distribution of v z for such ions0 The second 
term in Eq. (146 \ represents stimulated emission (treated as 
negative absorption), v' z is related to v z by the energy con- 
servation law m p Vz/2 + En,± + Tllu = m p v' z /2 + i?jv',±, an d 



3 This is essentially the partial cross section given by Eq. H33I I. 
except that the latter assumes v z = 0. To simplify notations, we 
suppress the subscript 'a' and the superscript 'ff' in this section. 



the integration is performed over those v z for which this law 
can be satisfied. In the case where J~n(v z ) = J- mp ,T(v z ) is 
the Maxwellian distribution (jTTJJ with T independent of TV, 
Eq. (146 1 can be written as 



H NN i(uj) = n N a NN ,(u) 



n N i j N > -N)0 p -Hw/T 



n N 



where a NN i(ui) = J Tn^z) a NN i(v z ,u;)dv z 
Eqs. (gSJ - (HZJ) yield 



LTE / n 

p (w) 



n v a{ui) (1 



-hu/T\ 



(47) 
In LTE, 

(48) 



where n p = ~^2 N tin and <t(u) is the average photoabsorption 
cross section of a proton: a(u>) = ~}2 N ffj ~}2 N , a NN i(cj). 

The power of spontaneous emission of unit volume into 
du dn is jui du ) dn, where jui is the emission coefficient. Some 
authors fe.g.. lZheleznvakovlll996l) call it emiss ivity (whereas 
other authors, e.g., iRvbicki fc Lightman|[l979l call emissivity 
the emission power per unit mass). It can be derived from 
the second Einstein relation (122[) and presented in the form 



/] 3lu,N'N, 
N'N 



3u),N'N 



hu/_ 

87T 3 C 2 



n N i I J r N >(v' z )a NN ,(v z ,iv)dv' z . 



(49) 



(50) 



In the case of Maxwell-Boltzmann distribution of v z and v' z , 
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The Coulomb logarithms for photoabsorption processes 
N — > N' strongly increase at ui fa (N' — N)ui cp . At these 
frequencies the weight of such transitions increases, which 
causes weak spikes (pse udoresonances) in the photo absorp- 
tion cross sections (see IPotekhin fc Chabrierl I2003T ). How- 
ever, at the p, T, and B values shown in Figs. [7] and [H] 
the resonant peaks of the Coulomb logarithms A*'^. 1 ^ at' 
(iV' > N) are not sufficiently high to make them larger 
than Aq'^q n . Therefore the free- free absorption/emission 
processes with N' — N give the main contribution even at 
uj w (N' — N)lj cp . Accordingly, these pseudoresonances are 
not very pronounced. They are not visible in the opacity 
curves in the upper panels of Figs. [7] and [5] because of the 
logarithmic scale, but the corresponding periodic decreases 
of the relative emissivity at multiples of uj cp are clearly seen 
in the lower panels. 
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Figure 8. The same as in Fig. UJ but for p = 100 g cm . In 
this case the relative emissivity for the weaker field (B = 10 14 G) 
equals to one because of the LTE. 



(51) 



using Eq. (|47[) we obtain 



Ju,N'N 



87T 3 C' 



I \ IN' -N)3„-hw/T 



In LTE, equations (gSJ , (g9j) and ([51]) reduce to the Kirchhoff 
law (for each polarization mode) 



jj LTE (uj)B U j(T)/2. 



(52) 



According to Eq. (|51[) . the ratio of the emission coefficient 
to its LTE value (often also called emissivity) is 



= E 



fi NN >(uj) , nui/T 



n N e (N-N')t3 p +hu,/T 

n N , 



(53) 



We shall call ratio (|53[) relative emissivity. 



6.2 Proton cyclotron line 

Figures [JJ and [8] show the opacities p(w)/p (upper panels) 
and relative emissivities j^/ii™ (lower panels) as functions 
of the photon energy for polarization a = +1, for two values 
of p and two values of B. At ui <C ijJ cp , the main contribu- 
tion to the absorption and emission of photons is given by 
the free- free processes preserving iV (N = N'). In this case 
~ ji TE - At higher uj, transitions N —> N' N give a 
noticeable contribution to the photoabsorption. In the ab- 
sence of statistical equilibrium, they result in a decrease of 
the relative emissivity (j^ < jt TE )- 



7 GENERALIZATION FOR OTHER IONS 

The formulae derived in the present paper for protons 
can be generalized for other nuclei with arbitrary A and 
Z. If they have spin |, then in the formulae for the 
rates of non-radiative collisions it is sufficient to replace 
the mass m p by mi = 0.9928 Am p , the number density 
rip by rii — n e /Z, the magnetic field parameter f3 p by 
0i ee huj ci /T = 73.9 (Z/ A) B 14 ,/T a , and to use for the De- 
bye screening wave number the Z-dependent formula = 
47r(l + Z)n e e 2 /T. Furthermore, in Eq. @ one should replace 
I NN ,(t/2) by I NN ,(t/2Z) and use Eq. (fX2Tjl . u' 2 = u 2 + 



2Z(N — N')m,/mi. Also in Eq. JT6]) one should substitute 
w(u±) and w x (u+, u_) from Eqs. (|A20|) (with Zi = Z 2 = Z) 
and (|A"24| . and to use u 2 = u 2 + Z(N - N' + N 2 - iVa ) from 
Eq. (HHJ. 

Generalization for ions with different spin is also 
straightforward, but more elaborate, because it requires to 
rewrite Eq. (|18p with allowance for different projections of 
spin on the magnetic field. 

A possible generalization of the free-free cross section 
a s for Z 7^ 1 is discussed at the end of Appendix [Bl 



8 SUMMARY 

We have derived the general expressions for the rates of tran- 
sitions between ion Landau levels caused by non-radiative 
and radiative electron-ion Coulomb collisions and non- 
radiative ion-ion collisions. We have also obtained (in Ap- 
pendix [B| the formulae for free-free photoabsorption cross 
sections in strong magnetic fields with allowance for the elec- 
tron and ion (proton) quantization, which are much simpler 
than the previously known ones. 

On the base of the calculated transition rates we solved 
the equation of statistical equilibrium for protons in a strong 
magnetic field. Considerable deviation from the Boltzmann 
distribution over the proton Landau levels occurs at densi- 
ties p < 0.1 B14 g cm" 3 . At higher densities (lower mag- 
netic fields) non-LTE effects are negligible. Conversely, at 
lower densities (higher B) the excited proton states become 
depleted because of radiative decay. 

Nevertheless, even with strongly depleted populations 
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of the excited Landau states, the emissivity of the fully ion- 
ized plasma is not much suppressed relative to its LTE value. 
This is because the main contribution in the photoabsorp- 
tion is given by transitions which do not change Landau 
number N. For such transitions, the relative population of 
other levels is unimportant. 

Although we have performed calculations only for the 
proton-electron plasma, generalization of our results to other 
ions is rather straightforward (Sect. [7}. 

All results in the present paper are obtained in the 
Born approximation. A truncation used in older papers to 
eliminate the divergence at small velocities, inherent to the 
Born approximation, is now replaced by a smooth correc- 
tion. The thermally averaged electron-proton non-radiative 
and radiative transition rates are not sensitive to this cor- 
rection. Proton-proton transition rates, however, are sensi- 
tive, therefore their more thorough examination beyond the 
Born approximation would be desirable. Fortunately, at the 
considered physical conditions the proton-proton rates are 
not dominant, thus we think that our main conclusions are 
sufficiently robust. 

Our general expressions for the proton/ion Landau 
level transitions derived in this paper will be useful for 
studying the possibility and the conditions of proton/ion 
cyclotron line formation in magnetar bursts. The radiation- 
dominated bubble formed during the magnetar outbursts 
may be considered as a hot (T > 10 keV), optically thick, 
rarefied medium embedded in a strong magnet i c fiel d 
ijThompson fe Duncan! 1 19951 : IWoods fe Thompson! l2005h . 
Within the bubble, vacuum polarization dominates the 
dielectric tensor and scattering dominates the opacity. The 
bubble may also contains appreciable amount of ions ripped 
out of the NS surface during the outbursts . A neutron star 
atmo sphere code, such as that dev e loped in iHo fe Lai! l|200ll . 
2003) and Ivan Adelsberg fe Tail (2006), 

can be adapted 

to study radiative transfer in the bubble. It may be that 
depending on the total bubble energy and the location of 
energy release in the bubble, the characteristics of bubble 
radiation (such as ion cyclotron line strength, line emission 
vs. absorption) are different. If so, the burst spectra can 
provide a useful diagnostics for the energy dissipation 
mechanisms of magnetar outbursts. 
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Al Scattering off a fixed Coulomb centre 

Al.l Derivation of the cross section for Z = 1 

The motion of a charged particle in a magnetic field can 
be described by different sets of wave functions, corre- 
spon ding to different choices o f the electroma gnetic gauge 
(e.g., lLandau fc Lifshitzlll976l ). IVentural l| 19731 ) studied the 
scattering by a fixed Coulomb potential using the axially 
symmetric gauge o f the vector potential A = x r. 
I Pavlov fc Yakovlevl (1 19761 ) found that the Landau gauge 
(A x , A y , A z ) — (0,—By,0) facilitates obtaining a simpler 
representation of the scattering rate. In the latter gauge, the 
'good' quantum numbers are the Landau number N and the 
components k x and k z of the wave vector of the particle. As- 
suming Z — 1, one can write the coordinate wave function 
as 

= (L x L z )-^ 2 e ik ^ +ik ^x N (y + k x al), (Al) 

where 

xAv) = J^(l% *«<VM (A2) 

and H N (0 = (-1) N d N /d£ N is a Hermite poly- 
nomial. The functions Xn(v) are ortho- normalized 
and have the following n ice property (|Klepikovl Il954l ; 
iKaminker fc Yakovlevl 1 1 98 ll) : 



r 

J — < 



Xn(v - 9*Om/2)x JV , (y + q x a 2 n /2)e qyV dy 



— T , In 2 r? I0\ pi(JV'-^)arctan( 9B / 9x ) 



(A3) 



Here q\ 



Q x + cfii, an d In'n is the Laguerre function 



|Sokolov fc Ternovfl986h . defined as follows: if N' - N > 0, 
then 



L N > N (X) 



^/2(N'-N)/2 r N'-N 



-^/2„(iV'-]V)/2 



X 
N 



(X) 



VN\N'\(- 



l\ (N — l)\ (N' — N + l)V 



(A4) 



otherwise In'n(x) = (—i) N ~ I NN r(x); L s N (x) is general- 
ized Laguerre polynomial. 

Relation (|A3[l allows one to reduce the matrix element 
of the transition \Nk x k z ) — > \N'k' x k' z ) to the form 



M 



2-kL x L, 



-iq H (fe a , + fc^)aJ 1 /2 + i(]V'-JV)arctan((j a / 9x ) 



xI N , N (q 2 _a 2 n /2)V q dq y , 



(A5) 



where q x = k' x — k x , q z = k' z — k z , and 



APPENDIX A: TRANSITION RATES FOR 
COULOMB SCATTERING 

Here we present the non-relativistic formulae for the colli- 
sion rates of charged particles in an arbitrary magnetic field, 
derived in the Born approximation. 



V q =e< 



dre lq 



47re 
q 2 + kl 



(A6) 



is the Fourier transform of the potential. In this paper we 
use the screened Coulomb potential for the electron-ion non- 
degenerate plasma, for which k s equals the inverse Debye 
screening length. 

Averaging the specific transition rate {2-k /h)\M \ 2 8{E' — 
E) (where E and E is the initial and final energy) over k x 
and summing it over k' x and k' z we obtain the transition rate 
per particle on the level iV with initial longitudinal velocity 
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v z = hkz/m, in a volume V = L x L y L z : 

u/ fix _ VzQjfW {Vz) 
W NN i = — 

*") Ly+J-oe J-oo h 



(2 

to rn e V ^— ^ «' 



dlfe' 



(A7) 



where to = h 3 /e 4 m e is the atomic unit of time, 0jvjv/(i>z) is 
an effective partial cross section, 



fix / \ 
W NN ,{U±) = 



dt. 



(A8a) 



u± = [(u ± It') + 



'^2, ,,2,1/2 



it = fc z |a m , it s = fc s a m , (A8b) 
(A8c) 



u = |fc^|a m = x/u 2 + 2A- 2 AT', 

and w^rjv' should be set equal to zero when u 2 <2{N'-N). 
It is easy to check that e quations fl ATI) . (IA8| are eq uivalent 
to equations (12), (13) of lPavlov fc Yakovlevl (|l976l ). 



A 1.2 Classical limit 

The fu nction nSW' (m± ) belongs to a class of integrals stud- 
ied bv IPotekhinl (Il996|). Accord i ng to his equation (B21), 
based on lKaminker fc Yakovlevl ((l98 1'), in the semiclassical 
limit 

W%, (u±) « (1/2) (ul/2 + AT + A') 

x[(\/7V 7 - \/iV) 2 +tt^/2]" 3/2 (A9) 

x[(^V 7 + VN) 2 + u±/2]- 3/2 (A>1, A'>1). 

Using this approximation, one can demonstrate that equa- 
tions (|A7p . (|A8|l provide the correct classical cross sec- 
tion in the limit B — > 0, where A can be replaced by 
{p±a m /Ti) 2 /2, p = m« = ftfc, p5_ = p 2 + p 2 ,. For exam- 
ple, consider a scattering event without screening (k s = 0), 
where the particle moves in the xz plane at the angle 9 to 
the z axis before scattering and 9' = 9 + a after scattering, 
6 < 7r/2, and 0' < 7r/2. Taking into account that p' — p and 
{PzPz -P±P±)/P 2 = cos a 
we obtain 



1 - 2 sin 2 a/2, from Eq. 



fix / \ 
W NN ,(U-) 



h 



cos 6 cos 9' 



8p 3 sin 4 a/2 



(A10) 



When A' ^> 1, a sum over A' can be replaced by an integral. 
Then, according to equations (IA7() and (|A10|) . the effective 
cross section for transitions into a range of Landau levels 
between A' and N' + dN' is 



da = 



4m m a 



dA = jj- w(U- 



TO Trie v z w to m e VTi 

TYe 4 m (1 — cos 9 cos 9') sin 9' 



2Vp 3 (1 — cos a) sin 4 a/2 



■<W. 



In the cylindrically symmetric case, when # - 
Eq. (|A1 1 f) becomes 



da 



4m 2 w 4 sin 4 a/2 



dfi Q 



(All) 
and 9' — » a, 



(A12) 



where dfic is a solid angle element. This is the Rutherford 
formula. 



A 1.3 Cross section and average transition rate for 
arbitrary Z 

For arbitrary charges of the scattered particle, Ze, and the 
Coulomb centre, Zoe, one should replace a m — > a m \Z\~ L / 2 
in equations (|A1[I - (|A5|I and K, — > Z ZV q in Eq. I|A5|) . Then 



fix / \ 47ra m m 2 /r^rv^ 
VzT m e ^-f 



,(u±) 



(A13) 



where w^j^i , w± , and it' are given by equations (|A8P with 
modified scaling of k z and fci, it = \kz\a m \Z\~ 1/2 and 
w' = | A;^. | a m | ^| 1//2 . Besides, the Z-dependence of the De- 
bye screening parameter u s = k B a m should be taken into 
account. 

If the velocities v z = hk z /mi have Maxwellian distri- 
bution (llUp . then from Eq. (|A13fl we obtain the probability 
for one particle in the state A in unit volume to make a 
transition to the state N' 



{VzO-NN>{v z )) = 



4v2tt m 

to m e 



\Z\at K 



where A„ 



30 du 



fix 

NN' i 



(A14) 



(A15) 



(3 = Ti\Z\eB/mcT. The function 9(u ) is the step function, 
equal to 1 when u 2 + 2(N — A') > and otherwise. 

A 1.4 Validity range and correction 



The integral in Eq. (|A15() diverges when N — N' . This be- 
haviour, which is well known for collision rates in the Born 
approximation in one dimension, means nothing but viola- 
tion of this approximation for low velocities of the colliding 
particles. 

A convenient parameter of the magnetic field strength 
at atomic scales is -y B = bj(a { ZZo) 2 = h 3 B/(m 2 c \Z\Z0e 3 ), 
where at is the fine-structure constant. Born approximation 
is valid at \k z \, \k' z \ ^> \ZZo\e 2 m/h 2 , that is at u, u 3> Jb 1 ^ 2 ■ 
The most important effect of going beyond Born approx- 
imation is the suppression of the amplitude of the lon- 
gitudinal part of the wave function (the exponential in 
Eq. |A1| ~) near the Coulomb centre. At \k z \ — > this ampli- 
tude becomes proportional to \/jfc7[. Quantitatively, when 
In jb 3> 1, the ratio of the square modulus of the am- 
plitude at \k z z\ 7 i3 1 ' /4 <C u <C 1 relative to its constant 
value at u > 1 equals C\k z \% 2 /(m\ ZZo\e 2 ) = Cu^/ys, 
wher e C = 2% / In 2 7s [1+0(1/ In 7b)] (jHasegawa fc Howard! 
1 196 il l. Therefore, at u -> or u -> 0, |M| 2 in Eq. (|A7) 
becomes proportional to « or 11', respectively, which com- 
pensates the diverging factor 1/u in Eq. (JATsJ . 

In order to eliminate the divergence of collision inte- 
grals in Bor n approximation, similar to Eq . (IA15II . previ- 



ous a uthors (jPavlov fc Panovl 1 19761 ; iKaminker fc Yakovlevl 
ll98lT ) introduced a cutoff at the lower limit of integration 
for N — N' . Instead of the cutoff, we introduce weight func- 
tion g(u)g(u'), with 

-1/2 



g(u) = (l + 7s u Y 



(A16) 



Under the conditions 75 > 1 and /?7s 3> 1, Born approx- 
imation is valid for most of the velocity values that sub- 
stantially contribute to the integral in Eq. (|A15[) . The latter 
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condition can be written asT> ZqZ 2 e 4 m* /Ti 2 , which is the 
usual condition of the applicability of Born approximation 
in the non-magnetic case. 

Apart from Born approximation, which consists in ne- 
glecting the influence of the Coulomb potential on the lon- 
gitudinal part of the wave function (the exponential in 
Eq. [Alj L we have also employed the adiabatic approxima- 
tion, which consists in neglecting perturbation of the trans- 
verse part of the wave function (x N in Eq. |Alj ). For the 
continuum wave functions, both approximations are always 
valid at z — > oo , but may become inaccurate at the dis- 
tances from the Coulomb centre comparable to the Bohr ra- 
dius. The adiabatic approximation remains sufficiently ac- 
curate at small z provided that the parameter 7s, intro- 
duced above, is large. Lowest-order perturbation corrections 
and exact solution to the continuum wave functions in a 
strong magnetic field beyond the adiabatic approximation 
have been discussed, e.g., by IPotekhin. Pavlov &: Ventura! 
l| 19971 ). 



A2 Scattering of two charged particles 

A2.1 Scattering of different particles 

Let us consider Coulomb scattering of two particles with 
charges Z^e (i = 1, 2). In the Landau gauge fSect. lATjl their 
wave functions are 

V'iv i ,fc Xii ,fc z . 1 (rO = (L x L z )~ 1/2 e lkx ' iXi+lkz - i!Si 

x\Z t \ 1/4 XN (\Z i \ 1 ^(y i + k x , i a 2 a /Z i )') , (A17) 



where Xn(v) ^ s gi ven by Eq. (|A2|I . The excitation energy 
of the two particles is E = heBc -1 QZi\/mi + \Z 2 \/m 2 ) + 
miVz : i/2 + m 2 v Z}2 /2, where mi are masses of the particles, 
and v z ,i their longitudinal velocities. 

Consider a transition in which the quantum numbers 
of the particles change from Ni,k X: i,k z ,i to Nl,k' xi , k' zi 
(i = 1,2). Wave functions (|A17[) depend on x and z only 
through the plane-wave exponential factor, which results 
in conservation of the x and z components of the total 
momentum in the matrix element of any potential which 
depends only on the relative position r 2 — r\ of the two 
particles: k xl + k x2 = k' xl + k' x2 , k zl + k z2 = k' zl + k' z2 . 
Since Xn does not depend on k Zt i, we may choose the 
reference frame comoving with the centre of mass in the z 
direction, so that k z2 = — k z \ = k z . The number of final 
states in dk' x l dk' x 2 dk' z is (2ir)~ Li L x L z dk' xl dk' x2 dk' z = 
(2-ir)~ 3 LlL z (ik xA dk' x ^dE' m t /h 2 \k' z \, where m„ = 
m\m 2 j(m\ + m 2 ) is the reduced mass. 

Using Fourier decomposition of the interaction poten- 
tial V{r) = (27r) -3 f dqe~ lq r V q (where r = r 2 - n), and 
assuming L x and L z to be large, we can perform the integra- 
tion over xi, q x , z, and q z in the matrix element M for the 
transition \Ni, N 2 , k x ,i, k Xi2 , k z ) — > \N{, N 2 , k' Xyl , k' Xt2 , k' z ) . 
This integration fixes q x = ki, x — k' lyX and q z — k' z — k z 
in V q . Furthermore, using Eq. (|A3[) , we can perform integra- 



tion over 2/1 and y 2 . Then we obtain 

' r dx 2 e i(q * +k *<2- k '*. 2):l 



M 



dq y V q 



XIn ^ \ 2\Z7\ 



2 2 

<7±a m 

V 2"2 \ 2\Z 2 \ 



x exp 



/ k x ^i + k x ,1 k Xi 2 + k x o 
ia m I — I q v 



2Z X 2Z 2 

X e' [( Jv l- Jv l) si S nZ i + ( JV 2- JV 2) si S nZ 2] arctan( Sy / 9!C ) (A18) 

By averaging \M\ 2 over k Xi2 and summing over k' x2 , we ar- 
rive at 



<|M| 2 ) 2 = J 



\Z 2 \L y 



2tv 



2iyL 2 L y L 2 



My \V q \ 2 lMM 



2 2 

q ± a m 



1 



j2 ( gjag, 
N * N * \2\Z 2 \ 



This expression does not depend on k Xt i, and therefore it 
does not require another averaging. Finally, summation of 
the specific transition rate {2-k /K)(\M\ 2 ) 2 S(E' — E) over 
k^. 1 and k' z gives the partial transition rate from Ni,N 2 
to N[, N 2 for two particles in volume V = L x L y L z 



W, 



N t N 2 ; N[N!, = V V z CT NlN2 . N ' iN ^(v z ) 



L z m 



h s \k 



ii y 



sign k'^ 



L dk k {m 2 } 2 



2tt 



47t m* 
to m e 



„j 1 

ZlZ 2 2 — 22 W N 1 ,N 2 ;N< 1 ,Ni('U±), (A19) 



V u 



where 



Jo 



At 



(t + ' 



t 



t 



\ 2\Z 2 \ 



(A20) 



JV i JV i \2\Zi\ 

Here u± — [(u ± it') 2 + u 2 ] 1 ^ 2 , u = |fc z |a m , u s — k s a m , and 

u 2 + ( Nl - N[) + {Ni n'J V2 . 

mi + m 2 mi + m 2 

(A21) 

The last equation represents the energy conservation 
law. The transition is energetically forbidden when the 
expression in square brackets is negative. Note that we 
have redefined u and u compared to the definition used in 
equations (|A13|I - (|A15|) : now Z % enter Eq. (|A21|l . Clearly, 

W N 1 ,N 2 ; N[ ,N' 2 i u ) ~ W N[,N 2 ; N lt N 2 ( u ) = ^JVi , ; JVj , N 2 ( u ) ■ 

In addition, if Z\ = Z 2 , then w Nl tN2 . jyj ,n 2 i u ) = 

W N 2 ,N 1 ; iV^.JV; ( u )- 



A2.2 Scattering of identical particles 

The derivation of the transition rates for identical particles 
can be patterned after Sect. IA2.11 but with initial and final 
wave functions in the form [tpN 1 ,k :cl ,k zl {ri)^N 2 ,k !c2 ,k z2 (r 2 )± 
^N 1 ,k ;cl ,k zl (r 2 ) ipN 2 ,k m2 ,k z2 {ri)] /-y/2, where ipN,k x ,k z (r) is 
given by Eq. (|A17|I . The resulting partial transition rate 
from Ni , N 2 to A^i , N 2 , averaged over initial and integrated 
over final k x values, is 



w 



^!,iV 2 ;JV(,]V^ - 2 (Wn^Nx, N[,N 2 ± ^i,JV 2 ;JV(,JV^)j (A22) 
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where the sign + (— ) refers to the states with even (odd) 
total spin, W Nl)Na . N ^ is given by Eq. (|A19[1 , 



ur x _ 871" m* |Z| 3 x . . 



(A23) 



and 



U N 1} N 2 ; N' 1 ,N 2 ( U -J u +) 



dt 



l0 (t + u 2 + /\Z\)(t + u 2 _/\Z\) 
xI N , Ni (t/2) I N , N2 (t/2) I N , N± (t/2) I N[N2 (t/2). (A24) 
The latter function satisfies symmetry relations 

W N 1 ,N 2 ; N[ ,N' 2 ( U - i U +) = 
W N 2 ,N 1 ; JVj ,JV^ ( W - ' U +) = 
W N[,N 2 ; N lt N 2 ( u_ ' u +)- 



W N 1 ,N 2 ;N! L ,N 2 ( U +, U -) 



APPENDIX B: CROSS SECTIONS OF 
FREE-FREE PHOTOABSORPTION 

In order to calculate a cross section of the free-free absorp- 
tion in a magnetic field, it is important to take into ac- 
count the the ion-el ectron centre of mass motion effects, even 
though m p 2> m e . IPotekhin fc Chabrierl (|2003T ) performed 
quantum calculations of this cross section and de monstrated 
a cor respondence to the classical dielectric tensor (|Ginzburd 
The result can be written as 



(T Q (w) = 



uj 2 vg(u) 



m e c (u! + aoj ce ) 2 (oj — acj cp ) 2 + uj^v^bj)' ^ ^ 
vt{u) =J2fnf» E ( B2 ) 

n,N n',N' 

D a (uj) — (l + a ^l^j u v {u) + (l — a ^L^j v e (u) + u^,{oj). 

(B3) 

Here a = 0, ±1 is the polarization index, t£ is the effective 
frequency of the electron-ion collisions for a given photon 
frequency uj, ft, and f n are the fractions of the protons 
and electrons in Landau states N and n, erf. n N . n i N i (o>) is 
a partial photoabsorption cross section for a transition in 
which the electron and proton change their Landau quan- 
tum numbers from n to n' and from N to N' , respectively. 
Finally, v p and v e in Eq. (|B3[) are effective damping fac- 
tors for protons and electrons, respectivel y, not related to 
the electron-proton collisions (for example, iGinzburgl (|l970T l 
considers collisions of electrons and protons with molecules) . 
The derivation of the damping factor (|B3|I from the complex 
dielectric tensor of the plasma assumes v £ <C uj ce , v a <C uj ce , 
and y p <C uj cv . 



IPotekhin fc Chabrierl i|2003h calculated 



n,N: n'N 



assuming LTE. Following their approach without LTE, how- 
ever retaining the Maxwell longitudinal distributions (|I0[) . 
we present the result as follows: 



ff,Q / S _ * 

l 'n,N,n',N'\ UJ l - n 



4 / 2ir n e e 4 , g 



3 V m e T huj 



(B4) 



where for a = 



/ 2 «(w' 2 ) 



3 /""dt. 0, u * 
2 Jo u> 6 

J2(u'±u) 2 w^ N . ]n , N ,(u ± ), (B5) 



and for a = ±1 



A 



n,N: n'N 



du 



-/3,u 2 /2 



9{u' 2 



Em„ ( uj cp \ 2 ea 

(B6) 



(m p + m e ) 2 

2 



Here /3* = hu*/T = heB /m*cT — /3 p m p /m*, and the ar- 
guments u± = [(it ± it') 2 + u 2 ] 1 / 2 , u = |fc z |a m , U s = fc s Q.m, 
and 



2 . 2m» , 2m* , 2w 
m H (AT — N ) H (n — n ) H 



1/2 



(B7) 



have the same meaning as in Appendix [S] The functions 
w(u±) are defined as 



..e,+l 



^n,s,n' -\-l,s' — 1 



-Vnv n -i, s+ i >n > >s :{p, m^2)J , (B8a) 

W n,iV;n',iV'( M ) = ^ ^ A P Vn,s,n> ,„> - 1 (P, uV2) 

-VlV + Tv n , s+ i, n >,s>(p,uV2)] 2 , (B8b) 

W n,N 1 ,n',N'( u ) = ^ p dp ^Vt? +T «n, s , n '+M'-l (p, «a/2) 

(p,uv/2)] [TiV 7 (p,W2) 

-VJV + l< 1+1 , n , a - (p, u\/2)] , (B8c) 

where s = N — n and s' — N' — n' are the relative proton- 
electron orbital quantum numbers, and v n>s>n i >s i (p, x) are 
the scaled Fourier trans forms of the effective p otent ials de- 
fined in Appendix B of IPotekhin fc Chabrierl (120031 ) PI Due 
to the symmetry properties of these potentials we have 



,(°) 



n,N: n'N' 



' N'( U ) = I PdpVn,s,n',s'{P;tJ-V2 



(B9) 



One can demonstrate that the integrals (IB8|) are sym- 
metric with respect to interchange of their indices n, n' or 
N,N', and moreover, all of them coincide with one another. 
It follows that Eq. ()B6|I simplifies to 



«ff,±l ,o / ,\ 3 f X du -f3,u 2 /2 , /2x 

K,N;n',N'(J 3 *> U /< J *)=2 J V 6 ( " ^ 

X [<]v ; n', N '(«+) +<lv;„',iV'("-)] - (BIO) 



(1) 



, (u) is any of the integrals (IB8 



where w. 

The integral (jB10|) diverges at u) — > u cp (N' - N) + 
io C e(n' — n), which is caused by the failure of Born approxi- 
mation at slow electron-proton relative motion. As discussed 
in Sect. IA1.4I this failure can be cured by introducing cor- 
rection factors (|A16|) under the integral. 



4 In equation (B7) oflPotekhin fc ChabrieU l|2003h the factor VIW 
(a typo) must be V s\s'\. 



The functions 10 j^. / N , (u) can be presented as 

1 P OO f\ot\ Ji 

= 2 J o (t + Jf 2)2 Z>,»®%'»(!)> ( B11 ) 

This alternative representation can be obtained by passing 
from the cylindrical to Landau gauge i n the derivation of the 
free- f ree cross section (Appendix B of IPotekhin fc Chabrierl 
2003), taking into account Eq. (1A3[) and recurrence relation 
L- n {x)-L' n _ 1 (x)=LiT 1 {x). 

We see that w^ N . , N , (u) coincides with w n ,N-. n ',tP ( u ) 
given by Eq. (%20J at Z 1 = Z 2 = 1, Thus, for a given 
initial (n,N) and final (n',N r ) Landau quantum numbers, 
the effective rates of the electron-proton radiative and non- 
radiative transitions are determined, for any B, T, and u), by 
one-dimensional integrals involving just two universal func- 
tions w }) i „,{u). These functions are smooth and mono- 

n,N ; n ,N' ' 

tonically decreasing. At u <C 1 they tend to constants, ex- 
cept for the following cases: (i) if n' = n and N' = N, 
then w ( V N;n , NI (u) ~ lnu and w^ N;n , N/ (u) ~ u -2 ; (ii) if 
ri = n and N' = N ± 1 (ai N' = N and n' = ra ± 1), then 
w n°jv-n' n'( u ) ~ n * m?i > wnere n * = max(iV, iV') (or n* = 
max(n,n'), respectively). At u ^> 1, iwijL -,/ N i( u ) u~ 4 . 

The approxima tion of infinitely massive ions 
ijPavlov fc Panovll 19761 ) is reproduced by setting tu cp — in 
Eq. (jBl| and replacing I%/ N (t) by J2n> — o n CO — 1 in 
Eq. (pT|) . 

Since Eq. (|B1[) is classical (although the factors and 
v a need quantum calculation), we may extend it to Z > 1. 
In this case the right-hand side of Eq. ()B4|I should be multi- 
plied by Z 2 , and the coefficient A n JV ;n',jv' will be different. 
By analogy with the general case of non-radiative Coulomb 
collisions, considered in Sect. lA2TT1 the latter difference con- 
sists in replacing of I% N > (t) by I 2 NN , (t/Z) in Eq. ||BT[) and 
multiplying (N — N') by Z in Eq. (|B7|) . This simple general- 
ization to the Z 7^ 1 case is possible only in the adiabatic and 
Born approximations, which we use in this paper. Beyond 
these approximations, the effects of centre-of-mass motion 
of an ion-electron system affect the initial and final wave 
functions in a non-trivial way. Continuum wave functions 
with allowance for these effects have bee n so far calculated 
only for Z = 1 l|Potekhin fc Pavlovlll997r i. 
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